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1. Introduction 

Diffractive production processes, being pure quantum mechanical by their origin, have always been a source 
of interesting information and sometimes even of confusion for all theoretical approaches to high energy 
scattering. It is enough to mention that due to the AGK cutting rules pj the total cross section of the 
diffractive production is equal to the contribution of the first shadowing correction (with opposite sign) to 
the total cross section. 

The goal of this paper is to discuss these processes in the framework of high density QCD |3|, f|, ||]. 
Since the most simple and elegant approach, that allows us to discuss both the total cross sections and 
the diffractive production on the same footing, is the AGK cutting rules, we choose the particular version 
of the theoretical approach to high density QCD: the BFKL Pomeron calculus [||, |J. This approach is 
based on the BFKL Pomeron || and the interactions between the BFKL Pomerons 0, || & [k|, which 
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are taken into account in the spirit of the Gribov Reggeon Calculus 11]. The BFKL Pomeron calculus 
can be formulated at least in two beautiful forms: as the functional integral || and as the evolution 
equation for the generating functional (see Refs.|l2], 13, |l4|, |1| [L6|]) . Both of them reproduce the mean 
field approximation for the total cross sections (see Refs. M, [L7], [l^]) as well as leading to the mean field 
equations for the diffractive amplitude (see Ref. |l9|]). In both of them we can include so called Pomeron 
(see Refs. [22, 23, 15, 24, 25, |26|]) loops which were not taken into account in the mean field approximation 



and which are very essential for the processes of the diffraction production especially in the region of large 
masses of the diffractive system. However, it is not clear at the moment the interrelation between the 
BFKL Pomeron calculus and the first attempt to suggest the equations for the diffractive processes in the 
framework of the approach to high density is more general than the mean field approach (see Ref. [£7| ) • 

In this paper we discuss the diffractive production in the simplest version of the BFKL Pomeron 
calculus, namely, in the toy model in which we neglect the fact that the QCD interaction can change the 
sizes of interacting dipoles. Such system of colourless dipoles reduces to the BFKL Pomeron calculus in 
zero transverse dimension. The analytical solution for the generating functional has been found in this 
model (see Ref. 



and an earlier attempt to solve this problem in Refs. Jl6|, 29, pQ|). Using this solution 



as well as the AGK cutting rules we can develop the theoretical approach to the diffractive processes. The 
first attempt to build such an approach has been made in Ref. [j36[ , however it was done in the limited 
kinematic region. 

The paper is organized as follows. In the next section we discuss the single diffractive production 
in both the mean field approximation and in the approximation that includes the Pomeron loops. The 
equations that govern this process, are written and the analytical solutions are found. However, we are 
able to solve the problem only in the restricted range of mass for the diffractive system (M), namely, 



1 



« Y = In (M 2 /m 2 ) < 



1 

a s 



hi 



(1.1) 



where m is a typical mass of soft interaction (say m = lGeV) and as is the running QCD coupling 
as = {N c /i{)as- In this region we can neglect the fluctuations in the diffractively produced cascade of 



partons due to the Pomeron loops contribution |3T], |33|, [3J, g5[ p6 |. It should be stressed that the Pomeron 
loops are taken into account for the "elastic" scattering with rapidity 



Y - Y = In (s/M 2 ) > —\n(^l) 
' a s \aij 



(1.2) 



In the third section we consider the double diffractive dissociation in the kinematic region of Eq. (1.1) 
for both produced masses (Mi and M<i) but at very high energy 

, / sm 2 \ 1 , /iV c 2 \ . , 
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The fourth section is devoted to calculations of the survival probability for a process with a large rapidity 
gap (LRG) of the following type 



hadron + hadron — * (1-4) 
hadron + { LRG } + {system of fixed mass ( di-jet, Higgs meson ,...)} + { LRG } + hadron 

which we can call a central diffractive process. The value of the survival probability of this process is 
a hot question since this reaction is one of the best for measuring Higgs meson at the LHC ( see Refs. 

@i m 0> ©in)- 



In section five we compare our estimates with other approaches that are on the market and in conclu- 
sions we summarize the results. 



Our approach in this paper is very close to one that has been developed in Ref. [36|. Therefore, it is 
worthwhile to stress the points on which we differ: 

1. In both papers the range of produced mass is restricted by Eq. (|1 . 1|) , and, therefore, we can neglect 
the fluctuations due to the Pomeron loops in produced system. It results in using the mean field 
approximations ('fan', tree Pomeron diagrams) for the produced system in both papers. However, we 
did not assume that the total energies restricted (Y — Yq <C In ( ^ ) ) as it is done in Ref. (36|] 



to justify the use of Mueller-Patel-Salam-Iancu approach [31, [3^, 33 1. For region, given by Eq. (|L 



we use the exact solution of the BFKL Pomeron calculus that has been found in Ref.|2 



In our approach we explore the Kovchegov-Levin equation |l9| (see also Refs. |2^, 21] where this 
equation was generalized in the mean field approximation) for diffractive production in the mean 
field approximation. Being equivalent to AGK cutting rules in the BFKL Pomeron calculus in zero 
transverse dimension, it has solid foundation in high density QCD (see Refs. [19, 27|) and allows 



us to develop a more general formalism in calculating the diffractive observables. In this formalism 
the solution to Kovchegov-Levin equation for Yq in the kinematic range of Eq. ([O]) plays role of the 
initial condition for the general equation for the scattering amplitude in the rapidity range of Y — Yq; 

3. Even in section 5.2 where we apply our approach to the Mueller-Patel-Salam-Iancu range of energies 
we obtain more general formulae since we do not assume that Y — Yq 3> Yq as in Ref. 



One of our motivations to discuss diffractive processes stems from the interesting paper of Y. Hatta, 



E. Iancu, C. Marquet, G. Soyez and D. N. Triantafyllopoulos [27| in which the diffractive dissociation has 
been considered in high density QCD. In this paper a generalization of Kovchegov-Levin equation |l9| is 
proposed. We wished to check and to confront these equations with the simplest model that we have at our 
disposal: the BFKL Pomeron calculus in zero transverse dimension. Indeed, we found two features of the 



diffractive processes in our approach that support the main ideas of Ref. [27]: (i) a kind of factorization 
between wee partons (Pomerons) produced by the incoming particle at rapidity Yq and the process of 
creation of the diffractive system by these partons (Pomerons); and (ii) an impressive difference between 
the mean field approximation and the exact answer. 
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2. Single diffractive production 



2.1 The mean field approximation. 



2.1.1 General approach. 




Our approach to diffractive production is based on the 
AGK cutting rules Q. These rules use the unitarity con- 
straints s-channel, namely, 



2ImA BFKL (s,b) = 2N 



BFKL 



(S,b) 



G 



BFKL, 



b) (2.1) 



where Im A BFKL (s, b) = N BFKL (s, b) denotes the imaginary 
part of the elastic scattering amplitude for dipole-dipole in- 
teraction at energy W = \fs and at the impact parameter 
b. It is normalized in a way that the total cross section is 
equal to a tot = 2 / d 2 b N BFKL (s,b). G BFKL is the contri- 
bution of all inelastic processes for the BFKL Pomeron and 
f d 2 b Gi n (s, b). Therefore, Eq. (2.1) gives us the struc- 



Figure 1: The diagram for the cut BFKL 



ture of the BFKL Pomeron exchange through the inelastic 

,> , n , c- , o i ■ processes and it can be formulated as the statement that the 

Pomeron that illustrates Lq. (2.1). ^ 

exchange of the BFKL Pomeron is related to the processes of multi-gluon production in a certain kinemat- 



ics (see Fig. ||). This equation was proven in Ref. ||. Using Eq. ( |2.1| ), we can express all processes and, 
in particular, the processes of diffractive production in terms of exchange and interactions of the BFKL 
Pomerons and/or the cut BFKL Pomerons (see Fig. |2|). 




Y 





Figure 2: Several examples of the Pomeron diagrams that contribute to the diffractive production of the bunch 
of gluon with ln(M 2 /m 2 ) = Y — Yq ; diagram a) is the first diagram for the process of the difractive production, 
while diagram b) gives an example of more complicated contribution to the diffractive production in the mean field 
approximation; diagram c) is the example of the diagrams that we neglected in the mean field approximation as well 
as in the approach of this paper. 

The AGK cutting rules establish the relation between different processes that stem from BFKL 
Pomeron diagrams. For example, the simple triple Pomeron diagrams in Fig. |3| leads to three inelastic 
processes: the diffractive production of the system with mass ln(M 2 /m 2 ) = Y — Yq (Fig. || -A); the multi- 
gluon production in the entire kinematic region of rapidity Y — with the same multiplicity of gluons 
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as in one Pomeron ( Fig. |3| -B); and the multi-gluon production in the region Y — but with the same 
multiplicity of gluons as in one Pomeron only in the rapidity window Y — Yq while for the rapidity Yq — 
the gluon multiplicity in two times larger than for one Pomeron( Fig. [3| -C). The AGK cutting rules Q say 
that the cross sections of these three processes are related as 

a A : °B ■ °c = 2 : -4 : 1 (2.2) 

At first sight the cross section of the process B is negative but it should be stressed that one Pomeron also 
contributes to the same process and the resulting contribution is always positive. 




Figure 3: The AGK cutting rules for the triple Pomeron diagram: the diffractive production of the system with 
mass ln(M 2 /m 2 ) = Y — Yo (Fig. || -A); the multi-gluon production in the entire kinematic region of rapidity Y — 
with the same multiplicity of gluons as in one Pomeron ( Fig. [| -B); and the multi-gluon production in the region 
Y — but with the same multiplicity of gluons as in one Pomeron only in the rapidity window Y — Yq while for the 
rapidity To — the gluon multiplicity in two times larger than for one Pomeron( Fig. ^ -C). The zigzag line denotes 
the Pomeron exchange while the zigzag line crossed y the dotted line stands for the cut Pomeron. 

Fig. ||, Fig. |2| and Fig. ||| as well as Eq. (|2.2D allow us to understand the equation for the single diffractive 
production in the mean field approximation that has been written by Kovchegov and Levin [18] and that 
has the following form (see Fig. || for graphical representation of this equation) 



dN D (x 01l b,Y,Y ) aC F 



dY 



7T 



2 



d 2 x 2 (2-3) 

J '02 x 12 



A^(x 02 ,b + ^x 12 ,Y,Yo) + iV D (x 12 ,b + ^xo2,T,To) - iV D (x i, b, Y, Y ) 

+ N D (x 02 ,b + Ixi 2) Y, Y ) iV D (x 12 , ix 02 , Y, Y ) (A) 

-4iV D (x 02 , b + ixi 2 , Y, Yo) iVo(x 12 , b + ix 02 , Y) (B) 

+2 AT (x 02 , b + ix 12 , Y) AT (x 12 , b + ix 02 , Y)\ (C) 
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A 



B 



C 




Figure 4: Different Pomeron cuts contributing to the cross section of diffractive dissociation which lead to different 
terms on the right hand side in Eq. ( |2.3| ) 

The physical meaning of Np is the cross section of the diffractive produced masses M large than 
ln(M/m 2 ) > Y-Y or in other words 



N D 



xoi,b,y,Y ) = / dY' N SD (x 01 ,h,Y,Y') (2.4) 
JYq 

where N SD (xQi,b, Y, Y') is the cross section of the diffractive production at energy Y with the produced 
mass equal to ln(M 2 /m 2 ) = Y — Y' . The initial condition for Eq. (|2.3|) follows from the definition of N D 
and it has the form with the initial condition given by 

7V£(x ± ,b,Y ) = N D ( XL ,b,Y = Y ,Y ) = 0; with - b ' Y ' Yo) \y=y = iV 2 (x ± , b, Y ). (2.5) 

where No is the solution of the Balitsky - Kovhegov equation [O] and Njj(x, b; Y; Yo) is the diffraction 
dissociation of the colourless dipole with size x± at impact parameter b into the system of gluons with the 
rapidity gap larger than Yo at energy Y. At first sight, Eq. ( |2.3| ) contradicts the AGK relations given by 
Eq. ( |2.2|) (see Fig. |^), but this contradiction can be easily resolved if we take into account coefficient 2 in 
Eq. Q (see Refs. @ for more details as well as for a proof based directly on the dipole approach to 
high energy scattering). 

Despite the simple structure of Eq. ( |2.3| ) which is only a little bit more complicated than the Balitsky 
- Kovhegov equation, as far as we know there exists no analytical solution to this equation and there is the 
only attempt to solve it numerically |43|. However, this equation has a simple solution in the toy model 



(see Refs. [19, 44, |45|] ) which we are going to discuss. 

2.1.2 The BFKL Pomeron calculus in zero transverse dimension. 

The mean field approximation looks extremely simple in the BFKL Pomeron calculus in zero transverse 



dimension (the toy model [12, 13, uM) . In this model Eq. fl2.3|) has a simple form 



^Wl^W = N d(u;Y,Y ) + N 2 D (u;Y,Y ) - AN D (u;Y,Y ) N (u;Y) + 2iV 2 (n;Y) (2.6) 
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Eq. ( p.6[) gives a typical example of the Liouville equation which has a general solution in the form 
No (u;Y) = Nq (Y + 4>(u)). Therefore u dependence stems only from the initial conditions of Eq. (|2.5|) 
and for Nq (it; Y) we have the Balitsky-Kovchegov equation flTj in a simple form Namely, 



dN (u;Y) 
dT(l -^2)Y 



N Q (u;Y) - N$(u;Y) 



ic iu,± ) - i'OlljJ = I o) = 1 - u = 7 



with the initial condition [12] 

NiJu^) = N (u;Y= Y/ : 

where 7 is the amplitude at low energy. Eq. ( |2.7D has an obvious solution 

iV lc (n,y ') exp(^) 



N (u;Y) 



(2.7) 



(2. 



(2.9) 



1 + N ic (uX) (exp(^) - 1) 

where y = T(l — > 2) K. and A/j c are the initial conditions for the scattering amplitude at Y = Yq. 
To solve Eq. fl2.6| ) we can introduce a new function (see Ref. |jl]] for details) , namely, 

N in (u;Y,Y ) = 2N (u;Y) - N D (u;Y,Y ) (2.10) 



for which we have the same equation as Eq. ( |2.7D but with the initial condition given by the following 
equation: 

N£( U] Yo) = N in (u;Y = Y ,Y ) = 2N ic (u,Y ) - N?{u,Y Q ) (2.11) 



Therefore the solution to Eq. ( 2.10 ) with the initial condition of Eq. ( |2.1l| ) has the form (see Eq. (|2.9| 

N£(u,Y ) exp(y-y ) 



N in (u;Y,Y ) 



1 + n£(u,y ) ( ew (y-y ) - 1) 

(2jVo(n,y ) - N*(u,Y )) ex P (y-y ) 
1 + (2N (u,Y ) - N*(u,Y )) (e W (y-y ) - 1) 



2.12) 



and 



N D (u;Y,Y ) = 2N (u;Y) - N in {u;Y,Y ) 



(2.13) 



Since N n (u;Y,Yo) satisfies the same equation as Nq (u;Y,Yo) (see Eq. (P?fl)) but with a different initial 
condition (Eq. ( pUp instead of Eq. (gp), one can see that the solution has the following formfl9|] 



N in (u;Y,Y ) 



N (u-y + Ay) with Ay 



In (2 + L^expQV 



(2.14) 



Using Eq. fl2.8|) for the initial condition of Nq (u, Y = 0) and Eq. ( |2.14| ) we obtain the following answers 
for the scattering amplitude, total inelastic contribution and for the diffractive production cross section 
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with mass in the interval \n(M 2 /m 2 ) < Y — Y : 



N (u;Y) 
N in (u;Y;Y ) 
N SD (u;Y;Y ) 

N D (u;Y;Y ) 



(1 - u) e y 
u + (l-n)e^ ' 

(l-u)e y (2 u + (1 -u)e y °) 
u 2 + (l-u)ey (2u + (l-u)e^)) ' 
u 2 (l-u) 2 e y + y ° 

(u 2 + 2u(i-u) e y + (i- u ) 2 ey+y°) 21 

2(l-u)e y (l-u)e y ( 2 u + (1 - u) e y ° 



u + (1 - u) e y 
in agreement with Refs. [19, 44, 45]. 



+ (1-u) e y (2u + (l-u)e^>) ' 



(2.15) 
(2.16) 
(2.17) 

(2.18) 



The advantage of our way of obtaining the solution is that it allows us to take into account an arbitrary 
initial conditions. We will use this in the more general case than the mean field approach. 



2.2 Diffractive production beyond the mean field approximation 

In the kinematic region of rather small produced mass (see Eq. (O)) but at very high energies (see Eq. (|1.2|)) 



we can neglect the correlations [33, 34, |36j in the diffractive produced system due to the Pomeron loops 
(see Fig. Q-c, for example). 



Y = Y 



mean field 
approximation 



Y=0 





mean field 
app roxitna t ion 



Y =0 



Figure 5: The general diagram for single diffraction 



production in the kinematic region given by Eq. (1.1) 
and Eq. (1.2) . The zigzag line crossed by the dot- 
ted line denotes the cut Pomeron while the same line 
without the dotted one stands for the Pomeron ex- 
change. 



Figure 6: The general diagram for double diffrac- 
tion production in the kinematic region, given by 
Eq. ( |l~T| ), for both masses \a{Mf/m 2 ) = Y - Y 1 
and ln(M|/m 2 ) = I2 and in the energy region , 
given by Eq. (Oh. The zigzag line crossed by the 



dotted line denotes the cut Pomeron while the same 
line without the dotted one stands for the Pomeron 
exchange. 
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The only diagrams that we need to consider are shown in Fig. || Our solution, given by Eq. (|2.12| ), sums 
these diagrams if we substitute in Eq. ( |2.12| ) the exact solution of the general equation for Ni n (u;Yq) = 
2 Nq (u; Yq) that includes the Pomeron loops. This equation has the form |]i~5| , |28|| 

dNo( U] y ) „ , ( dN (u;y ) ld 2 N (u;y ) 



du 

i)/r(i 



2). r(i 



du 2 
> 2) and T(2 



(2.19) 
1) describe the 



dy 

where we denote y = T(l —> 2) Y and k = T(2 
P — > 2P and 2P —* P transition, respectively. 

The initial and boundary conditions are 

Initial conditions: Nq (u; 3^o = 0) = 1 — it 
Boundary conditions: iVo (u = 0; y$) = 1 

Eq. ( 2.2C| ) stems from the fact that at low energies only the exchange of one Pomeron contributes to the 
scattering amplitude. The boundary condition at u = 1 says that the dipole at high energies does not 
interact with the target if this dipole interacts with zero amplitude at low energies ( Yq = 0). At u = the 
dipole interacts at low energy with the amplitude which is equal to the maximal value (N = 1) and any 
interaction between Pomerons cannot change this value of the amplitude. 



N (u = l;y ) 







(2.20) 
(2.21) 




. 8 



. 6 



. 4 



0.2 




2.5 



Figure 7: The amplitudes N (1 - u = 1/k, Y; Y ) and N in (1 - u = 1/k, Y\ Yq) as function of y - 34) = r(l — > 
2) • (Y — Yq) at different values of J^o 

The solution for the single diffractive production can be found within the following procedure: 

1. For the cross section of diffractive production of the mass M (Y = InM) the initial conditions for 
the solution of the master equation (see Eq. ( 2.19D ) we put at Y = Yo an d they have a form 

(I — u ) e yo 

N ^ = l + (l-u)(e»-l) ; (2 ' 22) 
Nm ^ Yo > ~ l + (i - 1) " [i + (l- u )(eM. - i)J ' (2 ' 23) 
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Figure 8: The total cross section of the diffraction dissociation in the range of mass (M) ln(M 2 /m 2 ) < Y 
No (1 — it = 1/k, Y; Yq) and the ratio of this cross section to the elastic amplitude ( the total cross section is equal 
to 2 N Q (1 - u = Y)) as function of y - y Q = T(l -> 2) • (Y - Y ) at different values of y 
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Figure 9: The cross section of the diffraction dissociation at in the range of 

mass y„,< + ^ (^o.i-^) = A^o determined by N SD (1 - u = 1/k; Y; Y 0ti ) 

(N D {l-u = l/K;Y;Y 0<i + ^) - N D (l-u = l/K-,Y-,Y 0ti - ^-)) /Ay with Ay a = 0.25 and the ratio 
N SD (1 - u = F; y 0ji ) /7V (1 - u = 1/k, Y") at different values of y 0ti . 



2. The next step is to find the exact solution (N$ xact (u, Y; Y ) and Nf* act (u, Y; In) ) to Eq. ( gig) 
with the initial conditions determined by Eq. ( 2.22j ) and Eq. ( |2.23 ). The general method for finding 



such solution is given in Ref. [28] and the only change that we need to make here is to include 
different initial condition. 



3. However, in Ref. [28| we use Kovchegov's formula [17, 14 1 for the amplitude, namely, 

N (l-u;Y) = 1 - Z(u;Y) 



(2.24) 
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where Z is the generating function. This formula assumes that n-dipoles can interact with one dipole 
(target) independently. It sums all diagrams of Fig. ^-b type. However 1 , if we would like to restrict 
ourselves by the enhanced diagrams of Fig. | we need to use a different relation between Z and N. 
Namely, 



Nq{1-u = 1/k;Y-Yq) 
N in (l-u= 1/«;Y;Y ) 



dZ(u;Y) 

du 
8Z m (u;Y) 



a 



u 



u=l I 
K 



\u=l 



(2.25) 
(2.26) 



In Fig. |7| the solution given by Eq. ( 2.25| ) and Eq. ( 2.26 ) is plotted as a function of Y — Yo at different 
values of Yq. 

4. Using Eq. (2.25) and Eq. ( |2.26| ) we can obtain the solution for the single diffraction based on 
Eq. ( 2.13| ), namely, 

N D (l-u= Y ; Y ) = 2 N (1 - u = l/«; Y; Y ) - N in {l-u = 1/k; Y; Y ) (2.27) 

and for the cross section of single diffraction in the system with mass Yq = In M we have 

dN D (l-u = l/K;Y;Y ) 



N SD (1 - u = Y; Y 



dY n 



(2.28) 



Our calculations for Njj (1 — u = Y; Yq) and N$d (1 — w = Y; Yq) are shown in Fig. ^ and Fig. ||. 
The most interesting observation is the fact that in a wide range of produced mass 3^o < 1 these cross 
sections only weakly depend on the value of the mass. 



3. Double diffractive production 

In this section we will discuss the double diffractive production, or the reaction of the following type 

{hadron (dipole)} + { hadron (dipole)} =^ (3-1) 
{ hadron(gluon) system with mass Mi} + {LRG} + { hadron(gluon) system with mass M2} 

where both masses (Mi and M2) satisfy the relation of Eq. ( |L2j ) while the total energy is so large that 
Eq. ( |1.3| ) holds. The diagrams that contribute to this process in the kinematic region, that we have 
mentioned, are shown in Fig. 0. For the double diffraction production in the kinematic region, given by 



Eq. (1.1), for both masses ln(M|ym ) = Y — Y\ and ln(M|/m ) = Y2 — and in the energy region, given 



by Eq. (1.3), we can neglect the Pomeron loops contributions in the region Y — Y\ and Y2 — (see Fig. g). 
This simplification leads to a simple procedure of calculating the cross section for the double diffractive 
dissociation, which consists of four steps: 



x We thank Edmond Iancu and Al Mueller for hot discussions on this subject. 



- 11 - 



1. As the initial conditions for the double diffractive processes we choose Eq. ( 2.22| ) and Eq. ( [2.23|) at 
Y = Y 2 . They give us N in (u, Y 2 ) and N (u, Y 2 ); 



2. Using the general method of Ref. [28], we find the generating function Z (u, Y\ — Y 2 ;Y 2 ) and 



Zi n {u,Y\ — Y 2 ;Y 2 ) as the solution of Eq. (2.19) with the initial conditions 



Z in (u, Y 1 -Y 2 = 0; Y 2 ) = 1 - N in (u, Y 2 ) ; 
Z(u,Y 1 -Y 2 = 0;Y 2 ) = 1 - N (u,Y 2 ) ; 



(3.2) 
(3.3) 



3. The third step is to notice that the amplitudes Nq (u, Y — Y\) and Ni n (u, Y — Y\) are also deter- 
mined by Eq. (|3|) and Eq. ( ^23|) at Y = Y - Yi 



4. The fourth and the last step is to find the final answer using the general formula (see Ref. [33 



13, [L4| ) that was re- written in the compact form in Ref. 35] 



n ( 7 = i/k,y -Yr,Y r ,Y 2 ) 

( d d 
exp — 7 



du\ du 2 
N in { 1 = 1/k,Y-Y 1 -Y 1 ;Y 2 
d d 

1 



exp — 7 



du± du 2 



N (Y - Yi,«i) Z(Yi;Y 2 ,u 2 j { \ Ul =l;u, i 
N in (Y - Y uUl ) Z in {Y r ,Y 2 , u 2 ) | \ U1 



(3.4) 
(3.5) 



1; U2=l 



We can simplify Eq. ( |3.4| ) and Eq. (3^) going to the Mellin transform for the generating functions 
Z in (Yi; Y 2 , u 2 ) and Z (Y 1 ;Y 2 ,u 2 ). Namely, 



ra+ico j 

Z(Y r ,Y 2 ,u) = / ^L^MiMz^.y^^ 



(3.6) 



Since d n No (Y — Y\,ui) jd 2 u\ is known, namely, 

d n N (Y -Y uUl ) 



i P y-yi x {^y-yi 



n\ e 



1) 



Eq. (|3.4D reads as 



N { 1 ,Y-Y 1 ;Y 1 ;Y 2 ) = 
a+ioo j p y-y\ 
' Z(Yi;Y 2l fi) 



(3.7) 



(3.8) 



2vri 



( e y-yi - i 



7T J \ 7 T 



exp 



7 T 



+ 1 



where T 



.y-yi 



l) and 7 = 1/k. Since (see Ref. |9| formula 3.381(3)) 



(3.9) 
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one can see that Eq. (|3.8| ) can be re-written in the form 

jv (7,y-ri;ri;r 2 ) x I ey _ yi j = (3.10) 

/ 1 \ f°° ( x \ r a+ico da 

- 7 Texp(_) ^ exp(-_) ^ (*; Y 2 , „) + 1 = 

-^expf^ xdxexp(-^\ Z (Y i; Y 2 , ln(l/x)) + 1 = 
1 \ Z 11 

Eq. ( |3.10| ) is convenient for calculation. In derivation of this equation we used that Z (Yi; Y 2 ,u = 1) = 1 
(see Eq. (ggg)). 

For iVj n (7, y — Yi; Yi; Y2) we obtain a reasonable approximation by using the same expression as 

Eq. fljU0|) but with T ^ 2T and Z — >■ Z in . Indeed, in Eq. (gjgg ) we can consider 1 — u ~ 1/k <C 1 and 

neglect the second term in Eq. fl2.23|) if (1 — u) exp(,y — ^o) *C 1 or 3^ — 3^o *C 1/(1 — u) ~ k. Recall that 

2 

5- 



7 Texp( — ) / _exp(-— -) Z(Y i; Y 2 ,u) + 1 



this interval for y — 3^o is much wider than the one given by Eq. (Fl), since k = l/a|. 



Finally, the inclusive double diffractive process, namely, the cross section defined as 

rY-Y! f Y 2 j 

Ng D = ^ dMMl/m 2 ) J q *MM$/m 2 )M?M$ (3.11) 

we calculate as 

N& D (u\Y,Yy,Y 2 ) = 2N (u,Y) - N in (u^Y^Y,) (3.12) 
For calculation of the double diffraction cross section in the two system of hadrons with fixed mass we need 

M 2 M 2 da _ d N® (u; Y, Yi , Y 2 ) 

MlM2 dMtdM% ~ dY,dY 2 (3 ' 13) 

4. Survival probability for the central diffractive production. 

In our approach the value, that we call survival probability, is defined as the following ratio: 

<IS 2 {> _ \A(F ig .m-a)\ 2 

<lSl> ~ \A(FigM-b)\ 2 { ' 

This ratio is the quantitative measure how the interactions of the 'wee' partons from different parton 
showers are suppressed the probability of the production of the central system, which is produced with 
rapidities close to Y/2 in the reaction of Eq. ( |1.4| ) type. 

To find the amplitude for the diagram of Fig. [lO|-a we need to take into account that the Higgs meson 



in Fig. 10 can be emitted from any of Pomerons at rapidity Y/2 in this diagram. We can take into account 
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Figure 10: The central diffractive Higgs meson pro- Figure 11: The diagrams that contribute in 

duction: diagrams Fig. |o|-a is a general case that Mueller-Patel-Salam-Iancu approximation, which 

takes into account the interaction of many parton describe the amplitude for the rapidities that satisfy 

showers, while the diagram of Fig. |Io|-b is the cen- the restriction Y = In (M 2 /m 2 ) < -J- In . 
tral production from one parton shower. 



this fact by the following procedure. First, we recall that the scattering amplitude can be written in the 
form \ 



14, 1, 



^ (-If d n Z(u,Y) ~ 



(4.2) 



n=l n=l 

Therefore, to calculate the Higgs production we need to take the following derivative 

dN( r ,Y) 

n=l 



7" 



(4.3) 



The next step is to solve our master equation (see Eq. ( 2.19 )) using Eq. ( [4 . 3[) as the initial condition. This 
solution (N (7; y ; y ) ) will allow us to find the survival probability. For calculating the survival probability 
we need to use the equation 

(N(r,Y-,m 2 



< \s z \ > 



7 exp (y) 



(4.4) 



In Fig. |T| our estimates for the value of the survival probability is plotted versus energy. We need 
to know the value of the P — > 2P vertex to assign the measured value of the total rapidity Y = ln(s/so) 
with so ~ 1 GeV 2 . Since this value is very important characteristic of the diffractive Higgs production we 
present the result of calculation in Table 1. If we take T(l — > 2) = Ap(0) where Ap(0) is the intercept 
of soft Pomeron (A P (0) « 0.1) the value of (|S 2 |) w 0.23. This value is in agreement with the estimates 
based on the two channel model for diffractive production 1 3*7 , [jjj. However, if Ap(0) ~ as = 0.25 the 
value of (|5 2 |) = 7 • 10^ 4 which is much smaller. It should be stressed that our calculations are one of the 
first that includes the enhanced diagrams (Pomeron loops) and can be used to evaluated the effect of this 
contribution on the value of the survival probability. One can see that the effect can be very large and the 
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< I s (JO I > 



0.8 
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r= 1/k = 1/16 
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J/ 



Figure 12: Survival probability (\S 2 ) versus y = T(l — > 2) Y at different values for the low energy amplitude 7. 
The solid line is our exact calculation (Eq. (4.4)) while the dotted line is the MPSI approximation (Eq. ( 5.13Q ) and 
the dashed-dotted line shows the calculation in the MFA (Eq. (|5.3|)) 



y 





0.25 


0.5 


0.75 


1. 


1.25 


1.5 


1.75 


2. 


2.25 


i\s 2 \) 


1 


0.96 


0.83 


0.75 


0.62 


0.751 


0.41 


0.31 


0.23 


0.16 


y 


2.5 


2.75 


3. 


3.25 


3.75 


4. 


4.25 


4.5 


4.75. 


5. 


{\s 2 \) 


0.11 


0.063 


0.045 


0.029 


0.017 


0.0077 


0.0054 


0.0034 


0.0018 


0.0007 



Table 1: The values of the survival probabilities versus y — T(l — > 2) Y at the value of the low energy amplitude 

7 = 



process of the diffractive Higgs production could be much less than it is expected based on the estimates 
in Refs. [37, 



39 . 



The estimates, discussed above, have one obvious defficientcy: they do not take into account the fact 
that the Higgs meson production cannot occur for energies less than s = M 2 where M is the mass if Higgs 
meson. In other language , the vertex for Higgs meson production has a size in rapidity while we neglected 
this effect in Eq. ( f4.4p . It is easy to take into account this replacing Eq. (4.4) by 



< \S 2 \ > 



with £ = In 



' M Higgs 
so 



(4.5) 
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y 


2.5 


2.75 


3. 


3.25 


3.75 


4. 


4.25 


4.5 


4.75. 


5. 


(\s 2 \) 


0.69 


0.39 


0.28 


0.18 


0.10 


0.048 


0.034 


0.021 


0.010 


0.0044 



Table 2: The values of the survival probabilities versus y = r(l — ► 2) y (the value of the low energy amplitude 
7 = 1/k.) Calculating (\S 2 \) Eq. was used with a s = 1/4 and M Higgs = lOOGeV. 

In the Table 2 we show values of < |S* 2 | > for 0:5 = 0.25, when calculations ware archived using 
Eq. ( |4.5| ). One can see that the difference between Eq. ( f4,4| ) and Eq. ( |4.5| ) (Table 1 and Table 2 correspond- 
ingly) is large and expected < |5* 2 | > at the LHC energy is equal to 4.4 • 10 -3 instead of (\S 2 \) = 7 • 10~ 4 
The conclusion that we want to make is that the enhanced Pomeron diagrams can be very important in 
obtaining the reliable estimates of the survival probability of the large rapidity gap at the LHC energy and 
can lead to additional suppression of the value of the cross section for diffractive Higgs production. 



5. Comparison with the approximate approaches. 



5.1 Our approach versus the mean field approximation. 



In the mean field approximation we can discuss the single diffractive production (see Eq. ( 2.15 ) - Eq. Q2.18D ) 



and the value of the survival probability. The direct application of the procedure described in Eq. ( |4.2[ ) - 
Eq. (fO|) for the MFA gives 

dN( T ,Y) = exp (|) 

dl (l + 7 (exp(^)-l)) 2 ' 

nUy-A = i^m . (5 . 2) 

V 2 J (1 + 7 (ex P (|) - l)) 2 + 7 ex P (|) (exp (|) - l) 



< \s 2 \ > = I 5-^ I ; ( 5 - 3 ) 

V(l + 7 (ex P (|) - l)) 2 + 7 ex P (|) (ex P (|) - l) J 

However, to compare the MFA with our calculations we need to go back from Eq. ( |2.25| ) and Eq. ( p.26| ) 
to Kovchegov's formulae 

N Q ( r ,Y) =1 - Z(l- T ,Y) ; (5.4) 
Ni n (7; Y) = 1 - Z in (l- T ,Y) ; (5.5) 



2 



The comparison with the MFA is given in Fig. [14] and Fig. 13 for the diffractive cross section and 
survival probability. As it is expected the MFA reproduces the exact result only for rather large values of 
the amplitude for dipole - target interaction at low energies. 
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Figure 13: Survival probability (|5 2 ) versus y = T(l — > 2) Y at different values for the low energy amplitude 7. 
The solid line is our calculation while the dotted line is the MFA approximation. 



5.2 Our approach versus Mueller - Patel - Salam - Iancu (MPSI) approximation. 



The MPSI approximation [31, 32, B3, 34] replaces the mean field one in the case if the dipole - target 



amplitude is small at low energies. The diagrams that contribute to the scattering amplitude in this 



approximation are shown in Fig. 11, they all have Pomeron loops but the size of any loop is larger than 



Y/2. The general formula of MPSI approximation can be written in the form 



N MPSI 



(7,Y) 



exp 



d d 



cxp 



r ye J 



du\ du2 
1 



N, 



MFA 




Y 

2 ,«1 



MFA 




— ,U2 



y 



r 0, 



7e 



y 



\U\=V ) U2 = l 

(5.6) 

r (0, x) is the incomplete gamma function (see formulae 8.350 - 8.359 in Ref. |4^]). We assume that 
e y/2 ^> 1 j n Eq, (^^) . In Eq. (|5.6j ) we substitute for Nq IFA the amplitude in the mean field approximation 
but we can use this equation using any of the amplitude calculated in Eq. ( 2.15|) - Eq. (|2.18 ). For example 
if one needs to calculate the single diffraction production Eq. (|5.6j) has the form 



N 



MPSI 
SD 



(l,Y,Y ) 



1 - < exp 
2 



d d 



du\ du2 



+ 1 



N 



MFA 
SD 



Y 



ui) N { 



MFA 




Y 



,U2 



1 + 27^ 



r 0, 



cxp 



where T x = T (Y,Y ) = 
In derivation of Eq. 



2 exp (y) - exp (y/2 + y ). 

) we use not Eq. ( |2.17| ) but simpler expression, namely 



|«l=l;tf2=l (5-7) 
(5.8) 



1 45] 



N, 



MFA 







(Y,Y ) 



2(1 



[2 e y 



1) 



(1 + (l-u) (2e3 



(5.9) 
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Figure 14: The ratio Nd/Nq as a function of y for different value of the amplitude of interaction for low energies 
7 = 1 — u. The solid line is our calculation while the dotted line is the MFA approximation. 



This simple formula does not take into account the elastic scattering of the dipole while Eq. ( [2.17 ) does 



Therefore, we find in Eq. (|5.7| ) the equation for the inelastic diffraction. Practically, such approach could be 
useful for the diffraction in hadron-hadron interaction, however, it is certainly is not good approximation 
for deep inelastic scattering where the elastic scattering also included in the experimental diffractive cross 
se3ction. 

Substituting in Eq. ( jO ) N^ IFA given by 



M FA = 2(l-u)e y ; w v r ? -— = (5.10) 

D v ' \\ + +(l-u) (e*-) 1 + (l-u)(2 e y - - 1) J y ' 

we obtain 

Ng PSI (r,Y,Y Q ) = (5.11) 
2|l- A; exp (-Vl T (V -V) - C f 1 - **p(-^A t4t 1 ( u - 



7^-^ \7 e ^/ V ' 7 e ^V Tj \ \7T1y 7T1 \ ' 7T] 
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Repeating the same calculation but for double diffraction using Eq. ([5.7D but with substitution Nq 
Nsd we obtain the following 



N^ SI ( 7 ,Y,Y U Y 2 ) 



—rT~ 



\{lT 2 ) 



1 + lT 2 



+ 1 



lT 2 V T 2 



x exp 



jT 2 



+ 



(lT 2 ) 2 lT 2 



where T 2 = T (Y, Yi,Y 2 ) = (2 exp (|) - exp(Ji)) (2 exp (£) - exp(^ 2 )). 

Finally, for the survival probability we obtain in MPSI approach the following expression 



< \S 2 \ > 



1 



1 + ~/e y 
H : — it-it exp 



1 



7 e' 



r o 



i 



7 e 



3-' 



(5.12) 



+ 4 



(5.13) 



In writing Eq. ( 5.13j ) we assume that 3> 1. 
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Figure 15: The comparison of the scattering amplitude for the dipolc- 
dipole scattering (No) calculated using Eq. ( [2.25| ). The x-axis is y — 
UBFKLOtgY — (4:ln2as/ '«) Y '. The solid line is the exact solution 
while the dotted one shows the MPSI result. 



For a comparison of the exact so- 
lution with the MPSI one we need to 
start from the dipole-dipole amplitude. 



In Ref. [28] we use Kovchegov's formula 



(see Eq. ( 2.24 )) for calculating the scat- 
tering amplitude. The key idea of this 
formula is that two or more dipoles at 
low energies interact with the target with- 
out correlations or in other words ac- 
cording to the Glauber approach. How- 
ever, for the case of the dipole-dipole 
scattering in the framework of pertur- 
bative QCD approach it was shown (see 
Refs. @, gj) that two or more dipoles 
do not interact with the target dipole at 
low energy. In this case we have to use 
Eq. fljTII ) and Eq. fl^ ) for calculating 



our amplitude. In Fig. one can see the result of these calculation. One can see from this figure that 
the MPSI approach describes quite well the high energy limit of the scattering amplitude but it underesti- 
mates the value of the amplitude at lower energies. This result is actually quite surprising since the MPSI 
approach was proven to be a good approximation in the limited range of energies, namely, 



1 

a s 



1 



— < Y < — In 



a s 



iV 2 



ai 



(5.14) 
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For our choice of k = 16 and as = 1/4 in Fig. |l5| ( considering u>bfkl = 4 In 2 as) Eq. ( [5.14 ) gives 



0.044 < y = {ojbfkl /k) Y < 0.2. One can see that in this range the MPSI approximation is not able to 
describe the exact solution. 



In Fig. 1£ we see that even for the single diffractive production the MPSI approximation describes the 



exact result quite well. 
6. Conclusions 

In this paper we obtain two main results. The first result of this paper is the set of formulae in the 
framework of the BFKL Pomeron calculus in zero transverse dimensions, that describe the diffractive 



dissociation processes in the limited range of produced mass given by Eq. (LI) - Eq. (|l,3|). The comparison 
with the approximation methods, that there exists on the market, shows that these methods are not able to 
describe the diffractive production in the entire kinematic region. However, the mean field approximation 
gives a very good description if the low energy amplitude of interaction with the target is rather large and 
if the two or more dipoles interact with the target without correlations according to the Glauber formulae. 

This great surprise stems from the fact that the MSPI approximation is able to describe the high 
energy limit of the amplitudes for elastic and diffractive scattering while theoretically this approach can 
be quaranteed only in the limited kinematic range given by Eq. ( 5.14]) . 



The surprise comes from two observations: first, the MSPI result does not lead to the decreasing cross 
section in spite of the fact that onle triple Pomeron interaction has been taken into account; and, second, 



for the total cross section the MSPI estimates are very close to the exact solution (see Fig. 15). 

Having in mind these two observations we believe that the MSPI approximation can be used as the 
simple first approximation instead of the MFA for the case of dipoloe-dipole scattering. However4, we 
would like to stress that for the value of the survival probability the MPSI approximation gives leads to 
higher probability of survival of the large rapidity gaps than the exact answer. 

The second result is the formula for the survival probability for the central diffractive (see Eq. (|1.4j) 
production of Higgs meson (di-jets and so on ). This formula is exact and it is valid for any value of energy. 
The estimates using this formula gives the value for the survival probability which is much smaller (in 10 
- 100 times) than the estimates of all approximation methods that have been used before. We hope that 
this result will stimulate a more detailed study of the survival probability in QCD and we consider the first 
papers on this subject (see Refs.pOj, 42] as only a beginning for such kinds of studies. 



In general, we hope that this paper, devoted to a model: the BFKL Pomeron calculus in zero transverse 
dimension, will be useful in searching the high energy asymptotic behaviour in the general case of QCD 
and/or in the case of the BFKL Pomeron calculus in two transverse dimensions. 

We would like to draw attention to the fact that in this paqper we sum mostly the so called emhanced 
diagrams (see Eq. ( 2.25| ) and Eq. ( 2.26| ) and Fig. [|,Fig. ^ and Fig. 10). These diagrams cannot give the 



correct answer since they did not reproduce the elastic cross section. The elastic cross section is equal to 
zero in this approach. However, accordingly to the general idea of the Pomeron calculus the sum of these 
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Figure 16: The comparison of the cross section of diffractive production with mass InM 2 > y^/ (ujbfkl a 2 s ) for 
the dipole-dipole scattering. The z-axis is y — y$ = lobfkl a 2 s (Y — Yq) = (4 ln2 as /a) (Y — Y ). The solid line is 
the exact solution while the dotted one shows the MFA result and dashed-dotted line shows the MPSI calculation. 



diagrams gives us the new resulting Green's function (Gp of the Pomeron and we should try to build a 
new theory of the interacting Pomeron based on this function. It is easy to see that at Y — > oo 

G P (Y,u) - > i (6.1) 

Therefore, the exchange of the resulting Pomeron leads (i) to the amplitude larger than unity but (ii) this 
amplitude is very close to the unitarity bound (Gp = 1). Therefore, we have a kind of black disc behaviour 
without the elastic cross section. We would lioke to recall that in the black disc limit a e i = \/ 2a tot- 

The fast decrease of the single diffractive production (see Fig. p|) shows that the effective triple Pomeron 
interaction vahishes at large Y. In other words, we can consider Gp given by Eq. ( |6.1[) as the exact solution 
of the multi interacting Pomeron system. On the other hand, the sum of enhanced diagrams gives the model 
for the two paricles irreducible diagrams in s-channel. To sum all diagrams we need to take into account 
also all two particle reducible diagrams. This sum is given by the eikonal (Glauber) formula, namely, 

Aasymp -» 1 - e~ Gp = 1 - exp {- - J < 1 (6.2) 

Therefore we obtain a grey disc limit. The eikonalization gives a simple solution for all observables calcu- 



lated here. However, Eq. ( p.2|) deserves a more careful consideration, as was noticed in Ref. [50], before 
being used. Eq. fl6.2|) gives the behaviour of the scattering amplitude at low energies but at u — ► (7 — ► 1). 
Therefore, solving the Pomeron calculus we reduce the problem of high energy scattering to the scattering 
at low energies while the last step: to find the scattering amplitude at u — * 1 remains open and, perhaps, 
eikonalization is one of possible solutions of this problem. 
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